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Abstract 

We construct sequences of axially symmetric multisphaleron solutions in SU(2) 
Yang-Mills-dilaton theory. The sequences are labelled by a winding number 
n > 1. For n = 1 the known sequence of spherically symmetric sphaleron solu- 
tions is obtained. The solutions within each sequence are labelled by the number 
of nodes k of the gauge field functions. The limiting solutions of the sequences 
correspond to abelian magnetic monopoles with n units of charge and energy 
E oc n. 
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1 Introduction 



Pure SU(2) Yang-Mills theory in 3 + 1 dimensions does not possess static, finite energy 
solutions. In contrast, SU(2) Yang-Mills-Higgs (YMH) theory possesses static, finite 
energy solutions, and so do SU(2) Einstein- Yang-Mills (EYM) and Yang-Mills-dilaton 
(YMD) theory. 

YMH theory with a triplet Higgs field contains a stable static, spherically symmetric 
solution, the 't Hooft-Polyakov monopole IQ , whereas YMH theory with a doublet Higgs 
field contains an unstable static, spherically symmetric solution |]^, ^, Q. This solution 
represents the electroweak sphaleron in the limit of vanishing Weinberg angle ^. 
For large Higgs boson masses the theory contains in addition a sequence of sphaleron 
solutions (without parity refiection symmetry) p, |^. 

EYM theory possesses a sequence of unstable static, spherically symmetric solutions 
for any finite value of the coupling constants p], Within the sequence the solutions 

are labelled by the number of nodes k of the gauge field function. When the dilaton 
field is coupled to the system, the solutions persist in the resulting Einstein- Yang- 
Mills-dilaton (EYMD) theory [|lll [T|, |13|, 0, |l6l. Decoupling gravity leads to 



YMD theory with the corresponding sequence of unstable static, spherically symmetric 
solutions 0, 

In YMH theory, beside the 't Hooft-Polyakov monopole, there exist multimonopoles 
with magnetic charge m = n/g, where n > 1 is the topological charge or winding 
number and g is the gauge coupling constant. The spherically symmetric 't Hooft- 
Polyakov monopole has n = 1. Following a pioneering numerical study |19||, axially 



symmetric multimonopoles have been obtained analytically in the Prasad-Sommerfield 
limit [2U|. In this limit the energy of the multimonopoles satisfies the Bogomol'nyi 
bound E = AT[n{(^)/g. 

Similarly, beside the electroweak sphaleron, which carries Chern-Simons number 
Ncs = 1/2, there exist axially symmetric multisphalerons in YMH theory with Chern- 
Simons number Ncs = ''^/2 PT|. So it is natural to ask, whether analogous axially 



symmetric solutions exist also in EYM or YMD theory. 

In this letter we construct multisphaleron solutions and their excitations in YMD 
theory. The appropriate axially symmetric ansatz for the multisphaleron solutions is 
analogous to the ansatz in YMH theory [l^, |22|, ^ . Like the YMH multisphalerons. 



the YMD multisphalerons are labelled by an integer n, which represents a winding 
number with respect to the azimuthal angle (p. While (p covers the full trigonometric 
circle once, the fields wind n times around. For n = 1, spherical symmetry and the 
known sequence of YMD sphaleron solutions are recovered. 

For each value of n we find a sequence of axially symmetric solutions, which can 
be labelled by the number of nodes k of the gauge field functions, analogous to the 
spherically symmetric case. For the limiting solutions, obtained for k oo, we give an 
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analytic expression. The energies of the hmiting solutions satisfy a Bogomol'nyi type 
relation E oc n, which represents an upper bound for the energies of the solutions of 
the n-th sequence. 

In section 2 we briefly review the lagrangian, discuss the ansatz and present the 
resulting energy functional. In section 3 we exhibit the multisphaleron solutions with 
n < 4 and k < 4. We discuss the limiting solutions in section 4 and present our 
conclusions in section 5. 



2 Axially symmetric ansatz 

Let us consider the lagrangian of YMD theory 

C = ^(9^$5^<f ) - ^e'^^'^TiiF^.F^n (1) 
with dilaton field $, SU(2) gauge field and field strength tensor 

F,. = d,V, - d,V, - ig[V,, K] , (2) 

and the couphng constants n and g. This theory possesses a sequence of static spheri- 
cally symmetric sphaleron solutions labelled by an integer /c, which counts the number 
of nodes of the gauge field function [|r^, |T^. 

To obtain static axially symmetric multisphaleron solutions, we choose the ansatz for 
the SU(2) gauge fields analogous to the case of multimonopoles |2^ and electroweak 
multisphalerons We therefore define a set of orthonormal vectors 

= (cosn0, sinn0, 0) , 
= (0,0,1), 

''^3^\4') = (sinn0, — cosn0, 0) . (3) 
and expand the gauge fields (V^ = V^t"-/2) as 

Vo^r-) = , V^rl = uf\c^)ut\<i^)w]{p, z) , (4) 
whereas the dilaton field satisfies 

$(f) = $(p,z). (5) 
Invariance under rotations about the 2;-axis and parity reflections leads to the conditions 
w\{p, z) = wKp, z) = wKp, z) = wKp, z) = wKp, z) = . (6) 
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The axially symmetric energy functional 

E = + Ev = j{e^ + e^'^^Ev) d(f)pdpdz 
contains the energy densities 



1 r 



and 



{dpwl + -{nwf + wl) - gwfwlY + {d,wl + -wl - gwlwlY 
P P 



+ {dpwl + -wl + gwfwiy + {d,wl + gw^wlf + {dpwl - d.wlf 



It is still invariant under gauge transformations generated by |19, 21 



(7) 



(8) 



(9) 



(10) 



where the 2-dimensional scalar doublet {wl, w^ — n/gp) transforms with angle 2r(p, z), 
while the 2-dimensional gauge field {w\, w'^ transforms inhomogeneously. We fix this 
gauge degree of freedom by choosing the gauge condition |^ 



dpwl + d,wl = 



'111 



Changing to spherical coordinates and extracting the trivial 6'-dependence (present 
also for spherically symmetric case n = 1) we specify the ansatz further pi] 



wl{r,e) 

wlir,e) 



— (l-Fi(r,9))cos9, wl(r,9) = - — (1 - F2(r,9)) sm9 



n 



(l-F3(r,^^))cos^^ , 

gr 



wl{r,i 



n 



— {l-FAr,9))sin9 . (12) 
gr 



With Fi{r,9) = F2{r,9) = F^{r,9) = Fi{r,9) = w{r), $(r,^) = y?(r) and n = 1 the 
spherically symmetric ansatz of ref. |jl8| is recovered. 

The above ansatz and gauge choice yield a set of coupled partial differential equa- 
tions for the functions Fi{r, 9) and <I>(r, 9). To obtain regular solutions with finite energy 
density with the imposed symmetries, we take as boundary conditions for the functions 
Fi{r,9) and $(r,^) 



r = 
r ^ oo 

9 = 
9 = n/2 



F,(r,^)|,=o = l, 
Fi{r,9)\r=oo = F{oo) 

deFi{r,9)\e=o = 0, 
deFi{r,9)\e=^/2 = 0, 



t = l,...,A, a,$(r,^)|,=o = 0, 

i = l,...,4, <l>(r,^)|,=oo = $(oo) 
z = l,..,4, de^r,9)\e=o = 0, 
^ = 1,...,4, de^r,9)\e=^/2 = 



(13) 
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(with the exception of F2{r,9) for n = 2, which has deF2{r,9)\0=o ^ [^) where 
F(oo) = ±1 and $(cx)) = 0. The boundary conditions for the gauge field functions at 
infinity imply, that the solutions are magnetically neutral. A finite value of the dilaton 
field at infinity can always be transformed to zero via — $(oo), r — T-e"****^""). 

The variational solutions -Fi(/3r, ff) and $(/3r, 6*), lead to E^^ = f3~^E^ + f3Ey. Since 



the energy functional is minimized for /3 = 1, the virial relation [IS 



E^ = Ev (14) 

also holds for general n. 

We now remove the dependence on the coupling constants k and g from the dif- 
ferential equations by changing to the dimensionless coordinate x = gr/n and the 
dimensionless dilaton function (p = The energy then scales with the factor l/{ng). 

The dilaton field satisfies asymptotically the relation 

lim x^09' = D , (15) 

where D is the dilaton charge. The energy is related to the dilaton charge via 

477 , , , 47r , , 

E = — \im xV = —D . 16 

Kg ^-^oo Kg 

3 Multisphaleron solutions 

We solve the equations numerically, subject to the boundary conditions eqs. (|T3|). To 
map spatial infinity to the finite value x = 1 , we employ the radial coordinate x = . 
The numerical calculations are based on the Newton-Raphson method. The equations 
are discretized on a non-equidistant grid in x and an equidistant grid in 6, where typical 
grids used have sizes 50 x 20 and 100 x 20 covering the integration region < x < 1 
and < 6 < 7c/2. The numerical error for the functions is estimated to be on the order 
of lO-l 

The energy density e, defined by 

E = — / e(x)x^dx sin 9 d9d(f) , (17) 
Kg J 

of the axially symmetric multisphaleron solutions and their excitations has a strong 
peak on the p-axis, while it is rather fiat along the z-axis. Keeping n fixed and varying 
k, we observe, that the ratio of the maximum energy density ^max to the central energy 
density e{x = 0) remains almost constant. Also, with n fixed increasing k, the location 
of the peak of the energy density approaches zero exponentially. On the other hand, 
with fixed k and increasing n the peak of the energy density moves outward. The 
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central energy density e(x = 0) as well as the maximum energy density ^max are shown 
in Table 1 for the sequences n = 1 — 4 with node numbers k = 1 — 4. Also shown in 
Table 1 is the energy E. 

In the following we exhibit as one example the multisphaleron solution for n = 3 
and k = 3. In Fig. 1 we show the energy density e. In Figs. 2a-d we show the gauge 
field functions Fi, which go from -Fj(O) = 1 to Fi{oo) = —1, passing three times zero in 
any direction. Finally, in Fig. 3 we show the dilaton function ip. 



4 Limiting solutions for k ^ oo 

The known sequences of solutions often tend to a simpler limiting solution. For the 
spherically symmetric YMD sequence with winding number n = 1, the limiting solution 
for — > cxo is given by \TE\ 

Woo = , V'oo = - In f 1 + - 

V X 

and describes an abelian magnetic monopole with unit charge, m = 1/g. The gauge 
field functions Wk approach the limiting function Woo = nonuniformly, because of the 
boundary conditions at the origin and at infinity, where Wk = ±1. 

For the axially symmetric sequences with tt, > 1 we observe a similar convergence 
with node number k. As for n = 1, the numerical analysis shows, that the gauge field 
functions (-Fj)fc tend to the constant value zero in an exponentially increasing region. 
For (Fj)oo = 0, the set of field equations reduces to a single ordinary differential equation 
for v?oo 

(xVj'-2e^--^ = 0. (19) 
This yields the spherically symmetric limiting solution 

(Fi)oo = 0, (^oo = - In (l + ^) , (20) 

corresponding to an abelian magnetic monopole with n units of charge. Thus the 
limiting solution of the sequence is charged, whereas all members of the sequence are 
magnetically neutral. This phenomenon is also observed in EYM and EYMD theory. 
(A detailed discussion of the convergence there is given in ref. ||16|| .) 

To demonstrate the convergence of the sequence of numerical solutions for ipk to the 
analytic solution ipoo, we show the functions cpk for k = 1 — 4 together with the limiting 
function ip^o in Fig. 4 for n = 3. The k-th function deviates from the limiting function 
only in an inner region, which decreases exponentially with k. The value <^fc(x = 0) 
decreases roughly linearly with k. 
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The energy of the hmiting solution of the sequence n is given by 



ng Jo 



2 , 47r 
X ax = — n 
Kg 



(21) 



i. e. the energy satisfies a BogomoPnyi type relation, E (x n. This relation is in 
agreement with eq. (0), since D = n for the limiting solution. This limiting value for 
the energy, eq. (pT]), represents an upper bound for each sequence, as observed from 
Table 1. The larger n, the slower is the convergence to the limiting solution. 

Further details on the solutions and the convergence properties of the sequences will 
be given elsewhere . 



5 Conclusions 

We have constructed sequences of axially symmetric multisphaleron solutions in YMD 
theory. The sequences are characterized by a winding number n, describing the winding 
of the fields in the azimuthal angle 0, while the solutions within each sequences are 
labelled by the node number k of the gauge field functions. For n = 1 the known 
spherically sjTiimetric sequence is obtained. 

The multisphalerons have a torus-like shape. The maximum of the energy density 
occurs on the p-axis. With fixed n and increasing k the maximum moves inward along 
the p-axis, whereas with fixed k and increasing n it moves outward. 

For fixed n and k —>■ oo each sequence approaches an analytically given limiting 
solution. The limiting solution has vanishing gauge field functions and corresponds to 
an abelian monopole with n units of magnetic charge. Because of the coupling to the 
dilaton, the energy of this limiting solution is finite. It satisfies a Bogomol'nyi type 
relation with energy E oc n. 



The spherically symmetric YMD sphalerons possess fermion zero modes |2J]. Since 
the electroweak sphaleron and multisphalerons have also fermion zero modes p5|, |21|, |26[ , 
we expect zero modes to be present also for the YMD multisphalerons. 

The spherically symmetric sequence of n = 1 YMD sphaleron solutions can be con- 
tinued in the presence of gravity, yielding a corresponding sequence of EYMD solutions, 
which still depends on a coupling constant |]TT], |T2|, |T^, |T^, |T^, |T^. By continuity, we 
conclude, that the axially symmetric n > 1 YMD multisphaleron solutions also exist in 
the presence of gravity, representing axially symmetric EYMD solutions. In the limit 
of vanishing dilaton coupling constant, they should reduce to axially symmetric EYM 
solutions, and thus axially symmetric generalizations of the Bartnik-McKinnon solu- 
tions p[. The outstanding question then is, whether there also exist the corresponding 
EYM and EYMD black hole solutions. 
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n 


ft 




6[X — Uj 


^max 


^p[x — UJ 


1 


1 


0.804 


12.48 


12.48 


-1.711 


1 


2 


0.966 


636.2 


636.2 


-3.374 


1 


3 


0.994 


25352. 


25352. 


-5.154 


1 


4 


0.999 


968441. 


968441. 


-6.968 


2 


1 


1.336 


0.5457 


1.734 


-1.399 


2 


2 


1.773 


7.256 


23.08 


-2.430 


2 


3 


1.927 


77.25 


251.0 


-3.524 


2 


4 


1.977 


840.2 


2674. 


-4.673 


3 


1 


1.800 


0.192 


0.92256 


-1.279 


3 


2 


2.482 


1.629 


7.552 


-2.110 


3 


3 


2.785 


10.79 


49.94 


-2.951 


3 


4 


2.913 


69.88 


322.4 


-3.837 


4 


1 


2.231 


0.0941 


0.6541 


-1.209 


4 


2 


3.137 


0.6388 


4.132 


-1.943 


4 


3 


3.588 


3.249 


20.52 


-2.649 


4 


4 


3.808 


15.92 


99.91 


-3.387 



Table 1 

The energy E (in units of the energy density at the origin e{x — 0), the maximum 
of the energy density Smax and the dilaton field at the origin (fi{x — 0) are shown for 
the solutions with node numbers k — 1 — A for the sequences n — 1 — A. 
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£(p,z) 71=3 k = 3 




Figure 1: The energy density e{p, z) (in units of ^) is shown as a function of the 
dimensionless coordinates p and z for the multisphaleron with winding number n = 3 
and node number A; = 3 
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F;(p,,zJ u=3 k = 3 




Figure 2a: The gauge field function Fi{pc, Zc) is shown as a function of the dimen- 
sionless coordinates pc and Zc for the multisphaleron with winding number n = 3 and 
node number A; = 3. The coordinates pc and Zc are defined by pc = ^^sin^ and 
Zc — cos 6*, respectively. The contourlines indicate the nodes of the function. 
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Dilaton Function ^k(^) n = 3 
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Figure 4: The dilaton functions (/9fc(x, for the multisphalerons with winding number 
n = 3 and node numbers k = 1 — 4 are shown as a function of the dimensionless 
coordinate x. The dashed, the dotted and the dash-dotted hues represent the angles 
9 — 9 — 7r/4 and 9 — 7r/2, respectively. Also shown is the limiting function i^oo{x) 
(solid hne). 
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